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A quantum master equation model for the interaction between a two-level system and whispering-gallery 
modes (WGMs) of a microdisk cavity is presented, with specific attention paid to current experiments involving 
a semiconductor quantum dot (QD) embedded in a fiber-coupled, AlGaAs microdisk cavity. In standard single 
mode cavity QED, three important rates characterize the system: the QD-cavity coupling rate g, the cavity decay 
rate K, and the QD dephasing rate Jj_ . A more accurate model of the microdisk cavity includes two additional 
features. The first is a second cavity mode that can couple to the QD, which for an ideal microdisk corresponds 
to a traveling wave WGM propagating counter to the first WGM. The second feature is a coupling between 
these two traveling wave WGMs, at a rate (3, due to backscattering caused by surface roughness that is present 
in fabricated devices. We consider the transmitted and reflected signals from the cavity for different parameter 
regimes of {g,$,K,y^}. A result of this analysis is that even in the presence of negligible roughness-induced 
backscattering, a strongly coupled QD mediates coupling between the traveling wave WGMs, resulting in an 
enhanced effective coherent coupling rate g = \/2go corresponding to that of a standing wave WGM with an 
electric field maximum at the position of the QD. In addition, analysis of the second-order correlation function 
of the reflected signal from the cavity indicates that regions of strong photon antibunching or bunching may 
be present depending upon the strength of coupling of the QD to each of the cavity modes. Such intensity 
correlation information will likely be valuable in interpreting experimental measurements of a strongly-coupled 
QD to a bi-modal WGM cavity. 

PACS numbers: 42.50.Pq, 42.60.Da 



I. INTRODUCTION 

Recent demonstrations of vacuum Rabi splitting in sys- 
tems consisting of a semiconductor microcavity and a sin- 
gle quantum dot (QD) Ol, 0, yj] represent an important mile- 
stone in investigations of cavity QED in solid-state materials. 
In these experiments, the microcavity-QD system is incoher- 
ently pumped with an excitation beam at an energy above the 
bandgap of both the QD and surrounding semiconductor ma- 
terial (usually GaAs or some form of its alloy AlGaAs). This 
pump light is absorbed and generates carriers in the GaAs sys- 
tem that can eventually (through phonon and carrier scatter- 
ing) fill the states of the QD; under weak enough pumping 
conditions, only the lowest energy bound exciton state of the 
QD is appreciably populated on average. Radiative recombi- 
nation of the exciton state and the resulting spontaneous emis- 
sion is then modified by the presence of a resonant microcav- 
ity. When the cavity is of small enough volume, the coupling 
(g) between the QD exciton and the cavity can be large, and if 
the cavity decay rate K and QD dephasing rate y± are smaller 
than g, the system is said to be strongly coupled [4], in that 
the QD exciton and cavity mode are no longer truly separate 
entities but are instead bound together. In the experiments de- 
scribed in Refs. HHHID], me evidence of this strong coupling 
has been presented in the form of spontaneous emission mea- 
surements from the QD-microcavity system, which display a 
double-peaked structure, rather than the single peak associ- 
ated with either the cavity mode or QD exciton alone. This 
vacuum Rabi splitting JUS] is one signature of the strong cou- 
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pling regime in cavity QED. 

Applications of strongly coupled QD-microcavity systems 
to areas such as nonlinear optics and quantum information sci- 
ence |0 Si, [BUI will also require an ability to effectively 
couple light into and out of the microcavity-QD device. That 
is, rather than measuring the spontaneous emission of the sys- 
tem alone, it is also important to have access to the cavity's 
optical response (transmission or reflection). This is true if, 
for example, one wants to examine the effect of a coupled 
QD-cavity system on the propagation of a subsequent beam 
through the cavity J3, E2L or if one wants to use the phase 
of the emerging transmitted signal within some type of logic 
gate II 1311 . Indeed, in most cavity QED experiments involv- 
ing an atom coupled to a Fabry-Perot cavity, it is the cav- 
ity's transmitted or reflected signal that is typically observed 

Q El EMI. 

Following demonstrations of coupling to s;7/ca-based cav- 
ities such as microspheres lfTsi [l9ll and microtoroids i20ll . we 
have recently shown that optical fiber taperslfl8l l2~lll are an ef- 
fective means to couple light into and out of wavelength-scale, 
semiconductor microcavities such as photonic crvstals j22ll23Tl 
and microdisks 1 24, 25]. In addition, we have shown that mi- 
crodisk cavities are extremely promising candidates for semi- 
conductor cavity QED experiments, with recent demonstra- 
tions of cavity quality factors (Qs) in excess of 10 5 |S!|26l] 
for devices with a mode volume (V e ff) of ~ 2 — 6(k/n) i . 
These Qvalues are significantly larger than those utilized in 
Refs. |EJ, [2, [H, and as a result, the devices that we con- 
sider are poised to operate well within the strong coupling 
regime, where multiple coherent interactions between the QD 
and photon can occur. It is envisioned that initial experi- 
ments in this fiber-coupled microcavity-QD system will ex- 
amine vacuum-Rabi splitting through measurements of the 
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FIG. 1: Illustration of the system under investigation. The 
microcavity-quantum-dot system is driven near resonance by cou- 
pling light into and out of it using an optical fiber taper waveguide, 
with a cavity-waveguide coupling rate K e (K e is a field amplitude de- 
cay rate). Imperfections in the microdisk cause a coupling of the 
clockwise and counterclockwise whispering-gallery modes, at a rate 
(3. These two whispering-gallery modes have a quantum-dot-cavity 
coupling rate go ar >d intrinsic cavity decay rate K, . The quantum dot, 
approximated as a two-level system, has a radiative decay rate Y|| and 
a total transverse decay rate y± • 



cently by Aoki, et al., in Ref. In this work our interest 
is to study this system in a parameter regime relevant to ex- 
periments involving the interaction of a single self-assembled 
semiconductor quantum dot with the microdisk WGMs in the 
presence of roughness-induced backscattering [33], and to de- 
termine the spectral response of the system for varying de- 
grees of quantum-dot-cavity coupling (go), backscattering (p), 
and modal loss (Kj). We examine how the phase and magni- 
tude of the backscattering parameter affect the coupling be- 
tween one or both cavity modes and the QD, and how the QD 
itself serves to couple the cavity modes together resulting in 
an enhanced coherent coupling rate over that of traveling wave 
WGMs. 

The organization of this paper is as follows: in section 
we review the simple classical coupled mode theory for modal 
coupling in microdisk cavities in absence of a QD. Sectionlllll 
presents the quantum mechanical analysis of this system in the 
presence of a QD. We review the quantum master equation for 
this system and look at semiclassical approximations for spe- 
cific choices of the backscattering parameter. In section [IV] 
we present the results of numerical solutions of the quantum 
master equation for parameters that are accessible in current 
experiments. Finally, the intensity correlations in the reflected 
cavity signal for various parameter regimes are studied in sec- 
tion |V] 

II. MODAL COUPLING OF TWO 
WHISPERING-GALLERY MODES DUE TO SURFACE 
SCATTERING 



transmission spectrum past the cavity; such measurements 
will be directly analogous to recent measurements of vacuum 
Rabi splitting from one-and-the-same atom in a Fabry-Perot 
cavity II 1711 . 

The goal of this paper is to provide a theoretical basis, 
accompanied by numerical simulations, for the experiments 
to be performed with single QDs in fiber-coupled microdisk 
cavities. Of particular concern is the proper treatment of 
the whispering-gallery modes (WGMs) in the cavities. More 
specifically, the WGMs have a degeneracy of two as modes 
with azimuthal number ±;« have the same frequency, but cir- 
culate around the disk in opposite directions. The WGMs are 
typically excited through an external waveguide, and for a 
nearly phase-matched system the forward propagating mode 
through the waveguide excites only the co-propagating mode 
in the resonator (the clockwise (CW) traveling wave WGM 
from here on out). Imperfections in the resonator will change 
this, as they cause backscattering that can couple the CW 
and counterclockwise (CCW) propagating modes (Fig. [TJ 
|ZZl HI IH HI. If the loss rates in the system (due to 
material absorption, scattering loss, etc.) are low enough, the 
backscattering can lead to coherent coupling of the CW and 
CCW modes, producing a pair of standing wave modes. A 
similar theoretical model focused on cooled alkali atoms cou- 
pled to dielectric whispering-gallery-mode microcavities has 
been presented by Rosenblit, et al., in Ref. |3jJ, and more re- 



The modal coupling between CW and CCW traveling wave 
modes in a whispering-gallery-mode microcavity has been ob- 
served experimentally and explaine d by many other authors, 
including those of Refs. [24j [27] [28j [29j [34j Here, we present 
a simple analysis of this coupling. This analysis is essentially 
an abridged version of that which appears in a recent paper by 
Borselli, et al., in Ref. |3fJ. 

Maxwell's wave equation for the vector electric field in a 
microdisk structure is 

V 2 E- A(0 (e + 5e)^ r = 0, (1) 

where po is the permeability of free space, £° is the dielec- 
tric function for the ideal (perfectly cylindrical) microdisk 
and 8e is the dielectric perturbation that is the source of 
mode coupling between the CW and CCW modes. Assum- 
ing a harmonic time dependence, the complex field modes 
of the ideal (8e = 0) microdisk structure can be written as 
Ej(r,f) = Ey (r)exp(i(Ojf), where j is an index label includ- 
ing the azimuthal number (;«), radial order (p), vertical or- 
der (v), and vertical parity (odd or even for a cylinder with 
mirror symmetry). In the microdisk structures of interest the 
vertical height of the dielectric cylinder is typically a half- 
wavelength in thickness, and only the lowest order vertical 
mode is well localized to the microdisk. In this case the verti- 
cal order and parity can be combined to define the fundamen- 
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tal vertically guided whispering-gallery-modes of the disk as 
transverse electric (TE), with antinode of the in-plane (p,$) 
electric field components at the center height of the disk, and 
transverse magnetic (TM), with antinode of the vertical (£) 
electric field component at the center height of the disk. In 
what follows we will continue to use the TE and TM des- 
ignation when discussing the WGMs, however, it should be 
noted that due to the radial guiding of the modes in the small 
microdisks of interest to this work the WGMs are far from ac- 
tually transverse electric or magnetic, and contain significant 
longitudinal field components in the azimuthal direction. 

Solutions to eq. (fl]i with 8e 7^ (i.e., modes of the perturbed 
structure) are written as a sum of the unperturbed mode basis 



E(r,0 



= e - mot 



L«;(0E°(r). 



(2) 



Plugging into eq. (Q]l, keeping only terms up to first order, 
and utilizing mode orthogonality, we arrive at a set of coupled 
mode equations 



da k 



— + iA(ti k a k (t) = z'J2 P./*:" j(t) 



(0o J5£(E»(r))*.E^(rMr 
2 Je°\E° k (r)\ 2 dr 



(3) 
(4) 



Reference 1 30] presents a functional form for p in situations 
involving small surface roughness perturbation. Under weak 
scattering conditions an assumption is made that only each 
pair (common radial order, etc.) of localized, degenerate CW 
and CCW WGMs with azimuthal mode number ±m are cou- 
pled by the disk perturbation 8e. The complex electric fields 
of the CW and CCW WGMs are simply relatedJH, and can 
be written in a cylindrical (p,(j),z) coordinate system as 



E cw( r ) = (E° p (p,z),iE;(p,z),E°(p,z))e im \ 
E° cw (r) = (E°Jp,z) -^(p, Z ),^ (p )Z )) e -^. 



(5) 



In the case of high-g resonant modes, with a small degree of 
loss per round-trip, the CW and CCW WGMs are to a good 
approximation complex conjugates of each other, which when 
combined with eq. (0 indicate that the WGMs can be written 
with transverse (p,z) electric field components real and longi- 
tudinal ((])) components imag inaryJH (i.e., E°,E9,E° all real 
functions). The coupled mode equations for these traveling 
wave modes then read as 



dacw 
dt 
flflccw 
dt 



= -;Acoacw(f) + '1Pl e '^ccw(f)7 
= -/AcoaccwO) + /|p|e~^acw(0> 



(6) 



with p = |p|e'^ given by (in a basis with the transverse electric 
field components of the WGMs real), 



P = 



0)0 J(/5 £e +' 2 '"^) ((E») 2 - (E°) 2 + (E° z ) 2 )pdpdz 
2 27tfe ((E°) 2 + (E°) 2 + (E°) 2 )p</p fl 'z 



(7) 

Equation © represents the time evolution of the two mode 
amplitudes (flcw^ccw) of an isolated system, without loss or 
coupling to an external waveguide. For the experiments con- 
sidered in our work, the waveguide coupler will be an opti- 
cal fiber taper through which light is traveling in the forward 
propagating mode. Light coupled from the fiber taper will 
selectively excite the clockwise WGM of the microdisk struc- 
ture due to phase-matching. Following the formalism of Ref. 
|36]fhis waveguide-resonator coupling can be included through 
the addition of a term ks + to eq. (O, where k is a waveg- 
uide coupling coefficient and |s + | 2 is the input power in the 
external waveguide (the squared magnitude of the mode am- 
plitudes, |a CVVjC ew| 2 , are normalized to stored optical energy in 
the cavity). Loss is introduced to the coupled mode equations 
by use of the phenomenological field amplitude decay rate 
Kt, taken to be the same for both the CW and CCW modes 
(though in general this does not have to be the case). This to- 
tal field decay rate is broken into a contribution from intrinsic 
microdisk absorption and scattering loss (K/) and a contribu- 
tion due to coupling to the external waveguide (K e ), so that 
K7- = K; + K e . Assuming lossless coupling and time reversal 
symmetry it can be shown 13611 that \k\ 2 = 2x e . The coupled 
mode equations then read: 



da 



cw 



dt 



(k t + /Aco)acw(f) + /|p|e^accw(0 + (iy/2K e )s+ 



d<2 j™ =-(K T + iA<i>)a C cw(t)+i\$\e ^a cw {t), 



(8) 



Here, the phase of the coupling coefficient was (arbitrarily) 
chosen to be purely imaginary, corresponding to a single-pass, 
waveguide transmission coefficient of +1 in the lossless cou- 
pler case ir36ll . These two coupled equations can be rewritten 
as uncoupled equations in terms of the variables asw.i and 
flsw,2, which represent the standing wave mode amplitudes 



flsw.i = — ^ («cw + e'^accw) 
flsw,2 = — j= (flew — e'^flccw) 



(9) 



As mentioned above, for an ideal microdisk the field distribu- 
tions of mode amplitudes acw and accw have an azimuthal 
spatial dependence of e ±,m, ^ so that the field distributions of 
«sw,i and «sw,2 correspond to (up to an overall phase fac- 
tor) standing waves \/2cos(m(j) — ^/2) and v / 2sin(m(j) — ^/2), 
respectively, with the azimuthal orientation of the standing 
waves being determined by the phase Z, of the backscattering 
parameter. Here, and in what follows, we take the origin of 
the azimuthal axis ((j) = 0) to be centered at the QD. 

The transmitted (Pj) and reflected (Pr) optical power in 
the external waveguide can be determined in either the basis 
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FIG. 2: Normalized transmitted (solid line) and reflected (dashed line) signal for standing wave whispering-gallery modes, determined through 
steady state solution of the coupled mode equations given in equation JgJ, (a) (3/kj- = 8, Kj- /k,- = 3 (b) P/kj- = 1, Kj- /k,- = 3, and (c) P/K7- = 1, 
K7- /k,- = 20. g; = 3 x 10 5 in all cases. 



of (acw^ccw) or (asw,i ,asw,2); because of phase-matching, 
coupling between the external waveguide and WGM resonator 
occur directly through (acw, fl ccw) and it is most natural 
to solve for these quantities in the traveling wave mode ba- 
sis. With the phase of the coupling coefficient chosen as 
purely imaginary the transmitted and reflected powers are 
P T = \s+ + (i\/2K7)acw| 2 and P R = \ (iy/2yCja C cw\ 2 , respec- 
tively. Steady state solutions for the normalized transmitted 
and reflected signals from the cavity for a number of different 
parameters are shown in Fig. [2] For p > Kt (Fig. EJa)), we 
see the formation of a distinct pair of resonances, located at 
co w (Do ± P- These dips correspond to standing wave reso- 
nances that result from a backscattering rate (p) that exceeds 
all other losses in the system (K7O so that coherent coupling 
between the CW and CCW modes can take place. As we see 
in Fig. |2jb)-(c), for P ~ Ky, the resonances begin to overlap 
and are no longer distinguishable. 

For cavity QED applications, one very important conse- 
quence of the distinction between traveling wave and standing 
wave modes is in the effective volume of the mode (V e ff), as 
the peak electric field strength per photon in the cavity scales 
as 1/ VVeff. In particular, we recall the definition of V e ff as: 



Veff = 



Je|E(r)| 2 
max[e|E(r)| 2 



(10) 



Standing wave WGMs have approximately half the volume of 
the traveling wave WGMs, so that the coupling rate g between 
a single quantum dot and a single photon in a standing wave 
cavity mode is expected to be V2 times that when the quan- 
tum dot is coupled to a traveling wave cavity mode. This of 
course assumes the single QD is positioned at an antinode of 
the standing wave mode; alternately, if it happens to be posi- 
tioned at a node, the coupling rate g will be zero. 

These arguments again rely upon having a physical system 
in which the backscattering coupling between CW and CCW 
modes is sufficiently strong compared to all other loss rates to 
allow for coherent modal coupling and formation of standing 
waves. They have also neglected the effects that an embedded 
QD may have, due to both an introduction of additional loss 
and mode coupling into the system. In the case of a strongly 
coupled QD we might expect that standing wave modes can be 



maintained provided that the modal coupling rate p exceeds 
not only Kj but also the QD spontaneous emission rate yii and 
non-radiative dephasing rate y p . To verify our physical intu- 
ition and understand the system in better detail, we consider 
a quantum master equation approach ll37ll to take into account 
the cavity-QD interaction. 



III. QUANTUM MASTER EQUATION MODEL 

We begin by considering the Hamiltonian for an empty mi- 
crodisk cavity (traveling wave WGM resonance frequency co f ) 
with field operators acw and accw and mode coupling param- 
eter p, written in a frame rotating at the driving frequency CO/ 
(and for h = 1): 



Ho = Aco c /fl cw^cw + Aco^/fl J, cw flccw — Pflcw^ccw 

— P*fl^CW^ cw + '(^CW - ^*^™)' 

where Aco c ./ = co f — CO/. As in the coupled-mode equations of 
the previous section, the CW propagating mode is driven by a 
classical intracavity field E = /\/2k p P/„, where K e is the cavity 
field decay rate into the waveguide and P,„ is the input power 
in the external waveguide. From this Hamiltonian, the classi- 
cal coupled-mode equations without dissipation can easily be 
derived through an application of Ehrenfest's theorem. 

Modeling the QD as a two-level system, we add the term 
Hi to the Hamiltonian: 



Hi =Aco Ci /d + d_ + 2g (tfcw°- ~ «cw6+) 
+ ^o(accw 6 - _ «ccw^+), 



(12) 



where Aco u / = co„ — CO/, co fl is the transition frequency of the ex- 
citon state of the QD, and go is the coherent coupling strength 
between the QD exciton state and the traveling wave WGMs. 
Note that go has been assumed real, and to have the same 
phase for both CW and CCW WGMs in eq. (TJ). This is 
consistent with a choice of the azimuthal origin lying at the 
location of the QD and for a QD dipole polarization trans- 
verse to the $ direction, where the electric field components 
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for both WGMs are equal and real at the position of the QD (a 
WGM basis can also be chosen in which this is true for dipole 
polarization parallel to <j)). For a QD located away from the 
azimuthal zero or with a mixed transverse and parallel dipole 
orientation, go will be complex, having a different phase for 
the CW and CCW modes. In general, care must be taken 
to calculate go and p consistently when studying interference 
effects between QD dipole scattering and roughness-induced 
scattering. 

The equation of motion for the system's density matrix p 
can be found from the equation: 



dp 1 , 



(13) 



where the term Lp = (Li + L2 + Lj,)p allows for the inclusion 
of decay through cavity loss (at a rate K7- = (O c /2Q), quantum 
dot spontaneous emission (at a rate Y||), and phase-destroying 
collisional processes (at a rate y p ), which are of particular im- 
portance for quantum dots, as unlike atoms, they are embed- 
ded in a semiconductor matrix where electron-phonon scat- 
tering is non-negligible. In the zero-temperature limit (appli- 
cable to the experiments under consideration as they will oc- 
cur at cryogenic temperatures), these loss terms are given by 
MM: 



Lip — K7-(2flcwPflcw — flcw^cwp — pflcw^cw) 

+ Kr(2a CC wPflccw — ^ccw^ccwP — Pflccw^ccw) 

(14) 



L 2 p = --M26_p6 + -d + d_p -pd+d_) 



i-3P = y(6,p6 z -p) 



(15) 



(16) 



From the master equation, we can numerically calculate the 
steady state density matrix p.„ and relevant operator expec- 
tation values such as (a^w^cw)™' which will then allow us 
to determine the transmission and reflection spectrum of the 
coupled QD-cavity system using formulas that are analogous 
to those used in the classical model of section lU These cal- 
culations are the subject of the following section. For now, 
however, we consider what intuition may be gained by further 
analytical study of the master equation. Using the standing 
wave mode operators, 



flsw,i = — ^= (,acw + e s accw ) 
flsw,2 = — (flew — e'^flccw) 



(17) 



and writing p = |p|e ! ^, we take operator expectation values to 
arrive at: 



dt 



(flsw.i) = -i(A(d c i- |P|)(dsw,i) 



-go -j= {<*-) - K nflsw,i) + -^= 



^(flsw,2> = -i(A(d c i + |P|)(flsw,2> 



1 — e'^ E 
-go — 7=-(g-) -^r(flsw,2> 



V2 



V2 



j t (a-) = -(jAG) fl/ +Yx)(a_) + i|((d z dsw,i)(l+^) 
+ <Msw,2)(l-e^) 



dt 



Or = - 



\/2go((a_fl SW , 1 t )(l + «*) + <d_ flsw , 2 t )(l 

- %/2g ((a + asw,i)(l + + (d+fl Sw .2> (1 - 
-7||(1 + 

(18) 

where we have used [a+,d_] = 6 Z and jj_ = y^/2 + y p . In 
the new standing wave mode basis both the empty-cavity fre- 
quencies and the QD-cavity coupling strengths are seen to 
be modified by the presence of strong backscattering. For 
the low-frequency mode (co c — |P|) corresponding to field 
operator flsw.i, the effective coupling strength is gsw.i = 
go(l + e*) /y/2, while for the high-frequency mode (co e + |P|) 
corresponding to field operator flsw,2, the effective coupling 
strength is gsw,2 = <?o(l — e'^) /V2. These coupling strengths 
are thus dependent on the phase t, of the backscattering param- 
eter p and can be as large as V2go or as small as zero. This 
result is consistent with what one would expect intuitively; 
the superposition of traveling wave modes results in a pair of 
standing wave modes whose peak field strength (per photon) 
is V2 times that of a traveling wave mode. The two standing 
wave modes are phase shifted from each other in the azimuthal 
direction by Jt/2, and as a result, if the QD is positioned in the 
antinode of one mode (^=0, so that gsw.i = V^go), it is within 
a node of the other mode (so that gsw,2 = 0), and vice versa 
for the situation when Z=K. Of course, for large cooperativ- 
ity, C, between the traveling wave cavity modes and the QD, 
the modes (flsw.i , flsw,2) ma y no longer be a good eigenbasis 
of the system. In order to gain some insight into such situa- 
tions before moving on to numerical quantum master equation 
simulations, we consider below steady state solutions to the 
semiclassical equations. 

The semiclassical equations of motion, derived from the 
above equations by assuming that expectation values of prod- 
ucts of operators equal the product of the expectation values, 
can be solved in steady state to yield information about the 
cavity response as a function of drive strength and detunings, 
and are useful for understanding the linear and nonlinear spec- 
troscopy of the systemj^, [TlJ, l39l l40l |4lll . In the case of a 
single cavity mode coupled to a two-level system this leads 
to the standard optical bistability state equation (OBSE). We 
consider two examples from the microdisk model described 
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above with two cavity modes, one in which the scattering 
by the QD and the roughness-induced backscattering are in 
phase and the coupling between the QD and cavity mode are 
described by the simple intuitive picture above, the other in 
which the two processes compete and the system response is 
more complicated. We begin with the simplest case in which 
^=0 (the \ = 71 case is identical except the roles of flsw.i and 
«sw,2 are swapped). Defining the parameters, 



and solving eq. (fT8T l with t, = in steady state we arrive at 
the following expressions relating the external drive (Y) to the 
internal state of the cavity (X + ,X-): 



Y±Y|| 
4gl 



8l 
2k t j± 



Y = 



X+ = 



(flsw,i)|^=o 



X_ = 



(flSW,2)|^=0 



(19) 



X+ = 



W+W+i V ^ ^MW^J (20) 



X- 



1 



Due to the common phase of the backscattering and the QD mode coupling in this case, the net effect of the backscattering on 
the system response is simply to shift the resonance frequencies of the empty-cavity modes. As expected the QD couples to one 
standing wave mode with a cooperativity twice that of a traveling wave mode, and is decoupled from the other. 

In the general case both standing wave WGMs couple to the QD and obtaining an equation analogous to the OBSE for an 
arbitrary t, is somewhat algebraically tedious. As a simple example in which both modes are coupled to the QD we consider 
\ = k/2, which yields in steady-state: 



X+ = 



l+!'AtO d /Kr + |fi|/K r 




K T 



4C 



X- = 



(21) 



In this case the backscattering and the QD mode coupling 
are out of phase and in competition. The resulting sys- 
tem response is governed by the detunings (Aco„/, Aco e ;) and 
the relative magnitude of the normalized roughness-induced 
backscattering, (|p|/Kr) 2 , and the cooperativity, C. 

Finally note that in order to connect to experiment an input- 
output expression between the incoming optical signal in the 
waveguide and the optical transmission (or reflection) past 
the cavity into our collection fiber is desired. In the formal- 
ism presented in section HI] the transmission and reflection are 
given in terms of the traveling wave mode amplitudes. These 
amplitudes can easily be recovered from even and odd parity 



superpositions of X + and X_ (cf. eq. (fT7b). 



IV. SOLUTIONS TO THE STEADY STATE QUANTUM 
MASTER EQUATION IN THE WEAK DRIVING REGIME 

The quantum master equation (QME) presented in the pre- 
vious section is solved numerically using the Quantum Optics 
Toolbox ll42ll43ll . We begin by considering steady state solu- 
tions, and calculate the transmitted and reflected optical sig- 
nals from the cavity. As a starting point, we eliminate the 
quantum dot from the problem by taking the coupling rate 



7 





FIG. 3: (a) Scanning electron microscope (SEM) image of a fabricated microdisk device. The disk thickness is /=255 nm and the sidewall 
angle is 9 = 26° from vertical. The measured average diameter for this device (i.e., the diameter at the center of the slab) is 2. 12 /m\. (b) Finite- 
element-calculated |E| 2 distribution for the TE p= i m= n WGM of a microdisk with a diameter of 2.12 fjm at the center of the slab. For this 
mode, X ~ 1265.41 nm, grad ~ 10 7 , and for a traveling wave mode, V e g ~ 5.6(A,/n) 3 . (c) Typical measured normalized optical transmission 
spectrum of the TE p= i g WGMs of a 2 fjm diameter microdisk similar to that in the SEM image of (a). 



go = 0. As expected, the resulting solutions (not displayed 
here) are identical to those obtained using the classical cou- 
pled mode equations and presented in Fig. [2] Having con- 
firmed that the QME solution is consistent with the classical 
solution in the empty cavity limit, we move on to study inter- 
actions with the quantum dot. To connect these simulations to 
ongoing experiments we choose physical parameters consis- 
tent with our fabricated devices 11251 12611 . In these experiments 
the microdisk cavity is 255 nm thick, and has a sidewall angle 
of 26° as shown in Figure[3ja). The modes of these structures 
(Fig- 0b)) can be numerically investigated through finite- 
element eigenfrequency calculations using the Comsol FEM- 
LAB software JlH 5J, |H], and information about the effec- 
tive modal volume V e g (as defined in eq. (fTOt ) and radiation- 
limited quality factor Q lad can be obtained. For the purposes 
of this work we focus on modes of transverse electric (TE) po- 
larization, where the electric field lies predominantly within 
the plane of the disk, and we consider first order radial modes 
(p = 1) in the 1200 nm wavelength band, the wavelength re- 
gion of the ground state exciton transition in our QDs. 

As discussed in ref. 026J, finite-element method simula- 
tions can be used to calculate V e ff as a function of the aver- 
age microdisk diameter D avg . From V e ff, we can estimate the 
QD-photon coupling strength. For a QD located at a posi- 
tion of maximum electric field energy density and with exci- 
ton dipole parallel to the local electric field of the cavity mode, 
go = d • Eph/ft is given by|0 H 

go=^\[M^> (22) 

2x sp y 27tn J> V e ff 

where x sp is the spontaneous emission lifetime of the QD exci- 
ton. Consistent with what has been measured experimentally 
for self-assembled InAs quantum dots l47ll . we take x sp = 1 ns. 
Figure [4] shows a plot of go versus disk size for traveling wave 
WGMs, and we see that go /2k can be as high as 16 GHz for 
the range of diameters we consider. As discussed in ref. l26ll . 
the WGMs are well confined (Q md > 10 s ) for all but the small- 
est diameter disks (D avg < 1.5 /jm). We have confirmed this 



in experiments !^ l48ll . with Q as high as 3.6 x 10 5 measured, 
so that cavity decay rates Kt/2% of approximately 1 GHz can 
reasonably be expected. Such devices exhibited doublet split- 
tings that are on the order of AX = 10-100 pm (see Fig. 0c)), 
corresponding to a backscattering rate |P|/2ti= 1-10 GHz. In 
practical devices then, the roughness-induced backscattering 
and the coherent QD-cavity mode coupling rates can be of 
similar magnitude, and we thus expect the QME simulation 
results to be particularly helpful in interpreting future experi- 
mental data. 

Unless otherwise specified, in all of the QME simulations 
to follow we consider the weak driving limit. In this limit the 
steady state response of the system behaves linearly, with the 
internal cavity photon number <C 1 and QD saturation effects 
negligible. For the QD and cavity parameters of the microdisk 
structures described below this corresponds to input powers of 
about 10 pW. 



ft 




41 , , , , 1 , , , 

1 fs£ J- 2.S 3 3.S 4 4.S 5 5.5 6 
a ve rage m ic rod i s k diameter tu ml 

FIG. 4: Coherent coupling rate go/1% (with T sp = 1 ns) for travel- 
ing wave TEp = i m whispering-gallery modes of the microdisk struc- 
ture described in Fig. [3] with varying disk diameter. Calculations 
were performed using a fully vectorial finite-element method, where 
for each microdisk diameter the azimuthal number of the TEp = i m 
WGM resonance was adjusted to place the resonance frequency near- 
est X = 1250 nm. 
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FIG. 5: Steady state QME solution for the normalized optical transmission (top curves) and reflection (bottom curves) spectra for a QD coupled 
to a microdisk cavity under weak driving and for three different QD detunings: (a) Aco flc = 0, (b) Aco ac = (3, and (c) Aco af = —(3. Cavity and 
QD parameters for these simulations are {go, (3, K7- , K e ,Y|| , J p }/2n = {6,9.6, 1.2,0.44,0.16,2.4} GHz, with the phase of the backscattering 
parameter set to i; = 0. In these plots the additional black dotted line plots correspond to an empty cavity (go = 0) and the red dashed line plots 
correspond to a QD with no non-radiative dephasing (y p /2n = GHz). 
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FIG. 6: Standing wave modes in a microdisk for 
different phases of (3 showing how the low and high 
frequency modes are positioned with respect to a 
fixed QD. (a) |3 > = 0), (b) |3 < (% = %), and 
(c)P = i|P| (| = 7i/2). 



A. |3>go > (kt.T-l) 

The first situation we study is one in which the backscatter- 
ing rate p exceeds the coupling rate go, which in turn exceeds 
the cavity and QD decay rates Kj and y_|_. We choose P/2tc = 
9.6 GHz = 0), with go/2% = 6 GHz, K T /2% = 1 .2 GHz (cor- 
responding to Q = 10 5 ), K e /2n = 0.44 GHz (corresponding to 
a transmission depth of 60% for the empty-cavity standing 
wave modes), and x sp = 1 ns (Y||/23t ~ 0.16 GHz). The unper- 
turbed cavity frequency (i.e., the resonance frequency of the 
traveling wave modes) is fixed at co 6 = 0, and three different 
QD-cavity detunings, Aco flC = co u — co c ={0, p, — p} are consid- 
ered. For each value of Aco 0( ., we calculate the steady state 
transmission and reflection spectra (as a function of probe 
laser frequency to cavity detuning, Aco/ £ = (0/ — co f ) from the 
cavity in three different limits: (i) go=0; here, there is no QD- 
cavity coupling, and the response should be that of an empty 
cavity, (ii) go/2n=6 GHz, y p /2n=0 GHz; here, we neglect 
all non-radiative dephasing, which becomes a good approx- 
imation as the temperature of the QD is cooled below 10 K, 
and (iii) go/27C=6 GHz, y p /2n=2A GHz; here, we allow for 



a significant amount of non-radiative dephasing, correspond- 
ing to a QD exciton linewidth of 10 /jeV, which is consistent 
with what has been observed experimentally at temperatures 
of around 10-20 K[|47ll. 

The results of the steady-state quantum master equation 
simulations are plotted in Fig. The interpretation of these 
results is as follows: as a result of the modal coupling due 
to backscattering, which has formed standing wave modes 
through a superposition of the initial traveling wave modes, 
only the lower frequency mode of the doublet has any spatial 
overlap with the QD (see Fig. |6]for location of the QD rela- 
tive to the two standing wave modes as a function of £), and 
thus, we should only expect the low frequency mode to exhibit 
any frequency shifts or splittings. In Fig. |2a), with the QD 
spectrally detuned equally from both empty-cavity standing 
wave modes, we see asymmetric vacuum Rabi splitting due to 
coupling of the QD to the low frequency mode at co c - — p. In 
Fig. |3b), with the QD now on resonance with the higher fre- 
quency mode, coupling still only occurs to the low frequency 
mode detuned in this case by 2p. Finally in Fig. |5l c X me 
QD is on resonance with the low frequency mode, and is also 
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FIG. 7: Steady state QME solution for the normalized optical transmission (top curves) and reflection (bottom curves) spectra of a QD coupled 
to a microdisk cavity under weak driving and for three different QD detunings: (a) Aav = 0, (b) Acfl ac = P, and (c) Aca ac = — p. These plots 
are calculated for identical parameters as in Fig. [5]with the exception that the phase of the backscattering parameter P has been changed from 
f;=0 to Jj=7C (P/2tc=-9.6 GHz). 




FIG. 8: Steady state QME solution for the normalized optical transmission (top curves) and reflection (bottom curves) spectra of a QD coupled 
to a microdisk cavity under weak driving and for three different QD detunings: (a.)Aa ac = 0, (b) Aa ac = P, and (c) Aco ac = — p. These plots 
are calculated for identical parameters as those in Fig. [5]wifh the exception that the phase of the backscattering parameter P has been changed 
from \ = to \ = jt/2 (0/271 = j'9.6 GHz). 



spatially aligned with it, so that we see the familiar symmet- 
ric vacuum Rabi splitting of this resonance. We note that the 
frequency splitting, Q.r, is in this case 2v / 2go rather than 2go; 
this is consistent with the mode volume of the standing wave 
modes being one half that of the traveling wave modes. For 
^ = K (Fig. [7]) the results are the mirror image of those in Fig. 
|3J where now the high frequency mode is spatially aligned 
with the QD and exhibits frequency shifts and vacuum Rabi 
splitting. 

Finally, we consider an intermediate backscattering phase 
£, = Jt/2. Here, we expect both modes to have an equal (but 
non-optimal) spatial alignment with the QD (Fig. [6f c)). The 
results, displayed in Fig. [8] show that this is indeed the case. 
In Fig. 13a), for example, we see a symmetric spectrum, con- 
sistent with both modes being equally spatially coupled to the 
QD and equally (and oppositely) spectrally detuned from it. 
In Fig. [Hfb)-(c), we see that the spectra are no longer symmet- 
ric, as the QD is on resonance with the high frequency mode 
in Fig. [Mb), and with the low frequency mode in Fig. Hfc). In 
each case we see Rabi splitting about the mode on resonance 



with the QD and only a small shift for the detuned mode. The 
Rabi splitting between the peaks is no longer at the maximum 
value of 2y/2go, but at a value closer to 2go due to the spatial 
misalignment of the QD with the empty-cavity standing wave 
modes. 

Before moving on to study different parameter regimes for 
{go, P,K,yx}, we examine the cavity's transmission spectrum 
as a function of the spectral detuning of the QD (Aco flC ). In 
practice |H 0> [U, QD-cavity detuning is often achieved by 
varying the sample temperature, which tunes at different rates 
the transition frequency of the QD (due to its temperature- 
dependent energy bandgap) and the cavity mode (due to its 
temperature-dependent refractive index). More recently, gas 
condensation on the sample surface H49il has been success- 
fully used to tune the cavity mode frequency of a surface- 
sensitive photonic crystal microcavity. Such a method hsa re- 
cently been shown to be effective for the microdisks studied 
here owing to the field localization at the top and bottom disk 
surfaces and at the disk periphery 115 Oil . In Fig. [9] we plot the 
cavity transmission minima as a function of A(O ac for the pa- 
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FIG. 9: Position of the resonance dips within the transmission spec- 
trum of a QD coupled to a microdisk cavity under weak driving, as 
a function of QD-cavity detuning Aco flc . {go, P)Kr,lQ.,Y||>Yp}/27l = 
{6, «9.6, 1.2,0.44,0.16,2.4} GHz, so that the QD is spatially coupled 
to both standing wave cavity modes. 



rameter set studied above in Fig. [8] where the QD is spatially 
coupled to both standing wave modes of the microdisk. When 
the QD is far detuned from the standing wave cavity modes, 
we see the response of an essentially uncoupled system, with 
transmission dips at the bare QD and cavity mode frequencies 
(±P). In the center of the plot, as the QD is tuned through the 
bare-cavity resonances, a pair of anti-crossings are evident as 
the QD couples to each of the standing wave modes of the 
microdisk. 



B. go > P> (K r ,Y±) 

Here we switch regimes to one in which the QD-cavity cou- 
pling rate dominates all other rates in the system, including 
the backscattering rate p. In particular, we choose go/27i=12 
GHz, with (3/27t=4.8 GHz, Kr/27C=1.2 GHz (K e /2;t=0.44 
GHz), and T sp =l ns (j\\/2% — 0.16 GHz). The qualitative 
behavior that we expect to see is similar to that of the pre- 
vious section as both go and p represent coherent processes, 
so that their relative values are not as important as their val- 
ues in comparison to the dissipative rates in the system. This 
is seen in Fig. [TOta), where the QD is spectrally located at 
— p, so that it is resonant with the low frequency mode of the 
standing wave doublet. Predictably, the interaction with the 
QD causes this resonance to split, with a splitting Qn=2\/2go. 
The higher frequency mode remains unaffected, as the choice 
of Z,=0 causes it to be spatially misaligned with the QD. 



C. K r > go > P > n 

Now, we take the cavity loss rate Kr/2;r.=9.6 GHz to ex- 
ceed both g /27l=6 GHz and P/2;r.= 1.2 GHz. In addition, 
K e /2% = 3.5 GHz, Yii /27t=0. 16 GHz, and J p /2n=0 or 0.7 GHz, 



so that Kt > K e > Yx (g°°d cavity limit). In the absence of 
a QD we expect to see a single transmission dip rather than 
a doublet for Kr ^> p. This is confirmed in simulation by 
the black dotted line in Fig. [ToT b). With the addition of a 
QD, taken to be resonant with the center frequency of the sin- 
gle cavity transmission dip, we expect to see this single dip 
split into two, with the dips not being completely resolved 
due to decay of the cavity mode (K7- > g). This is confirmed 
in Fig. [ToTb), where the splitting £l R /2% = 14.8 GHz lies be- 
tween the expected splitting for a purely traveling wave cavity 
mode (Q.R=2go) and the expected splitting for a purely stand- 
ing wave cavity mode (fi«=2\/2go), and lies closer to the for- 
mer due to the large degree to which Kr exceeds p. 



D. Y|| > SO > P > k f 

Here, the roles of Kr and Y|| are swapped in comparison to 
the previous subsection, so that Y||/27t=9.6 GHz is the domi- 
nant dissipative rate, exceeding each of {go, p,K7-,K f ,}/27t = 
{6,1.2,0.6,0.22} GHz (bad cavity limit). Unlike our pre- 
vious example, in absence of a QD we do expect to see a 
pair of standing wave modes form, as p > Kr. This is con- 
firmed in Fig. [TOt c) (black dashed line). Now, we introduce a 
QD that is spectrally aligned with the low frequency mode 
at —p. Because QD decay is so large in this case we ex- 
pect that the standing wave character of the modes is going 
to largely be erased when coupled to the QD. To confirm this 
intuition, we examine the calculated transmission spectrum in 
Fig. [TOlc). The low frequency mode does indeed split, but the 
splitting Q.r/2k = 14.4 GHz is less than the expected splitting 
of 2\/2go for standing wave modes, and lies much closer to 
the 2go splitting for traveling wave modes. The situation thus 
mimics that of the previous example, although in this case the 
relatively weak transmission contrast of the QD-coupled res- 
onances is a result of operation in the bad cavity limit. 



E. go > K r > P > Y± 

Finally, we consider a scenario in which the QD-cavity cou- 
pling go/271 — 12 GHz is the dominant rate in the system, but 
where cavity decay Kr/27C=6 GHz exceeds the backscatter- 
ing rate p/27C=1.2 GHz. In absence of a QD we see a sin- 
gle transmission resonance dip (Fig. [Told)) as Kr > P- If a 
QD is now spectrally aligned to the center of this transmis- 
sion dip (Aco ac =0) three resonances appear within the trans- 
mission spectrum of Fig. [TOl d). This should be contrasted 
with the transmission spectrum of Fig. [ToTb) in which only 
two resonant transmission dips were present. The central res- 
onance dip of Fig. [Tol d) is at a detuned frequency of 1.2 GHz 
(= P/2JI), and corresponds to the frequency of one of the two 
standing wave modes that can form through an appropriate 
combination of the traveling wave modes. As this mode is 
spatially misaligned from the QD for t, = 0, we do not expect 
its frequency to have shifted due to interaction with the QD. 
The other two transmission resonances correspond to the split- 
ting of the low frequency mode from its empty-cavity position 
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FIG. 10: Steady state QME solution to the normalized optical transmission (top curves) and reflection (bottom curves) spectra of a QD coupled 
to a microdisk cavity with % = 0, and for: (a) go > P > K7- > Y± (Ad) m . = -(3/27t, {go,P,K r ,K c ,Y||,y p }/2n = {12,4.8,1.2,0.44,0.16,2.4} 
GHz), (b) k t > go > P > 7± (Aco flc = 0, {go, p, K T ,K e ,y h y p }/2n = {6,1.2,9.6,3.5,0.16,0.7} GHz), (c) Y|| > g > P > K r (Aa ac = 
-P/27T, {g ,P,Kr,K e ,Y||,Yp}/27t = {6,1.2,0.6,0.22,9.4,0} GHz), and (d) g > k t > p > Y± (Aaw = 0, { g0 ,$,K T ,Ke,y\\,yp}/2n = 
{12,1.2,6,2.2,0.16,0.7} GHz). In these plots the additional black dotted line plots correspond to an empty cavity (go = 0) and the red 
dashed line plots correspond to a QD with no non-radiative dephasing (y p /2n = GHz). 
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FIG. 1 1 : Steady state QME solution for the normalized optical trans- 
mission (top curves) and reflection (bottom curves) spectra for a 
QD coupled to a microdisk cavity under weak driving and with the 
roughness-induced backscattering rate (P) zero, go > Kj > Y± > 
P (Aco flc = 0, {g ,$,KT,Ke,Y\\,y p }/2n = {6,0,1.2,0.44,0.16,0.7} 
GHz). The additional black dotted line plots correspond to an empty 
cavity (go = 0) and the red dashed lines plots correspond to a QD 
with no non-radiative dephasing (y p /2% = GHz). 



at -p/27i = -1.2 GHz. The splitting of £l«/27i;=33.6 GHz is 
very close to the value of 2v / 2go expected for interaction with 



a standing wave mode. 



The basic result that the above example demonstrates is that 
the QD can effectively serve as a means to couple the traveling 
wave microdisk modes, even in instances where the backseat - 
ter parameter is small relative to other rates in the system. As 
a final illustration of this, we consider the situation where the 
backscatter parameter is zero. In Fig. QT| the empty-cavity 
single transmission resonance separates into three resonance 
dips, one at the original zero detuning and the other two split 
by The interpretation of this result is that the QD has 

effectively served to couple the two counter-propagating trav- 
eling wave modes, creating a pair of standing wave resonant 
modes, one which is decoupled and has an electric field node 
at the position of the QD, and the other which is strongly cou- 
pled to the QD at a field antinode. In this case, and in other 
strong coupling cases where go is the dominant system rate, 
the QD serves to set the position of the effective standing wave 
cavity modes (as opposed to backscattering phase c^) thus en- 
suring azimuthal alignment of the QD with a field antinode 
of one of the standing wave modes. Note that in the example 
of Sec. |IV A| (Fig. [5]) in which (3 > g , it is the phase of P 
(t) which determines the position of the standing wave field 
antinodes with respect to the QD. 
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FIG. 12: Normalized second order auto-correlation function (solid blue line) g C c W ,ccw(t) f° r the counterclockwise propagating cavity mode 
for the parameters {go,(3,K7-,K e ,Y||,Yy,}/23t = {6,9.6,1.2,0.44,0.16,2.4} GHz, Aco ne = — p, and whose transmission/reflection spectrum, 

originally shown in Fig. 0c), is re-displayed here in the plot insets, (a) CO/ = — P — goV2, (b) CO; = — P + goV2, and (c) CO/ = p. The additional 
red dashed line plots correspond to a QD with no non-radiative dephasing (Yp/27t = GHz). 



V. INTENSITY CORRELATION FUNCTION 
CALCULATIONS 

Of additional interest is the behavior of the (normalized) in- 
tensity correlation function g( 2 \x), whose value can indicate 
nonclassical characteristics of the cavity field [51], and is thus 
of essential importance in the characterization of QD-cavity 
based devices such as single photon sources 

Gamma. 

Furthermore, intensity correlations of the cavity field (and 
other higher-order correlations) are sensitive to the energy lev- 
els of multi-photon states of the system, and thus provide fur- 
ther information about the system beyond the weak driving 
limit studied above. This is particularly important in the case 
of WGM cavities, in which the presence of a double-peaked 
spectrum typically associated with Rabi-splitting cannot by 
itself be regarded as evidence for strong coupling. Here we 
analyze the intensity correlations of a coupled QD-cavity sys- 
tem, the cavity containing a pair of nearly-degenerate WGM 
modes as in the analysis of the previous sections. 

A general definition for any (stationar y) tw o-time intensity 
correlation function in our system is 11551 15611 : 



[a\t)tf(t + x)b(t + x)a{t)) 
P(t)a(t))(fc(t+x)b(t+x)) 



where d and b are the field annihilation operators for modes 
a and b, which can be the cavity traveling wave modes (la- 
beled cw/ccw) or standing wave modes (labeled sw\/sw2). 
Here, it is assumed that steady-state has been reached (i.e., 

(2) 

t — ► °°), so g a ^(t) is the stationary two-time correlation func- 
tion, and is a function of the time delay X only. We calculate 
g„L(x) by applying the quantum regression theorem 13711 and 
numerically integrating the quantum master equation (eqns. 

(fTTTi-(fT6]l) li42l l43ll . In what follows, we initially focus on 

(2) 

calculating ghcw,ccw(t), the two-time intensity auto-correlation 
function for the counterclockwise WGM field operator, flccw- 
Due to phase-matching, the reflected signal from the cavity is 
proportional to accw, allowing sucn intensity correlations to 
be measured in practice. 



We begin by considering the set of parameters studied in 
the steady-state transmission and reflection spectrum of Fig. 
[2c), where t, = so that p is purely real and positive, and 
where Aco flC — — p so that the QD is tuned to resonance with 

the empty cavity lower frequency standing wave mode which 

(2) 

it is spatially aligned with. g C cwxcw{x) is calculated in three 
instances, with each case corresponding to a probe field fre- 
quency CO/ tuned onto resonance with one of the three reso- 
nance peaks in the coupled cavity-QD reflection spectrum of 
Fig. |3c). The results are shown in Fig. [12] For probe frequen- 
cies CO/ = — p±go\/2 (Fig. [T2l' a)-(b)). photon antibunching 
and sub-Poissonian statistics are predicted. This antibunching 
is a result of the anharmonicity of the Jaynes-Cummings sys- 
tem; once the system absorbs a photon at — p±go\/2. absorp- 
tion of a second photon at the same frequency is not resonant 
with the higher excited state of the system 01211 . The degree of 
antibunching is a function of the specific system parameters 

f21 

chosen, and gccw,ccw(0) approaches zero more closely as go 
further exceeds the rates Kt and y± . In this case, the difference 

in g^ccw,ccw{x) for probe frequencies CO/ = — pigov^ is a re- 
sult of the asymmetry in the spectrum of the system due to the 
presence of the nominally uncoupled high-frequency standing 
wave mode (asymmetry in the probe frequency detuning from 
the CCW traveling wave mode's natural frequency also plays 

^) a role here, and has a persistent effect upon ghlv,ccw{x) even 
for increasing mode-coupling (p, go) and decreasing dephas- 
ing (K7-, Yx))- For a probe frequency resonant with the third 
reflection peak at CO/ = p (Fig. Etc)), g cc Lccw (x) is essen- 
tially unity for all times and the reflected light from the cavity 
is nearly Poissonian due to the spatial misalignment, and re- 
sulting de-coupling, of the QD from the high-frequency cavity 
mode light field. 

We next examine the parameter set explored in Fig. [3a), 
where t, = n/2, so that both standing wave cavity modes are 
spatially coupled to the QD. In addition, Aco oc = 0, so that the 
modes are equally and oppositely detuned from the QD. Once 

(2) 

again, we calculate g ccw ,ccw (x) for three cases, with each case 
corresponding to CO/ on resonance with one of the three peaks 
in the reflection spectrum of Fig. HJa). The results, shown in 
Fig. [13] indicate mild antibunching for CO//27T, = ±12.8 GHz 
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the parameters {go, PjK^Ke, Y|| , Yp}/2t = {6, (9.6, 1.2,0.44,0.16,2.4} GHz, Aco flc = 0, and whose transmission/reflection spectrum, originally 
shown in Fig. [8]a), is re-displayed here in the plot insets, (a) C0//2JI = —12.8 GHz, (b) (Hi/2% = GHz, and (c) C0//27t = 12.8 GHz. The 
additional red dashed line plots correspond to a QD with no non-radiative dephasing (y p /2n = GHz). 
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FIG. 14: Normalized second order correlation function (solid blue line) gcch>,ccw(t) f° r the counterclockwise propagating cavity mode for 
the parameters {go, P, K7-, K e , yn , Yp}/27t = {6,0, 1.2,0.44,0.16,0.7} GHz, Aco ae = 0, and whose transmission/reflection spectrum, originally 

shown in Fig.QT| is re-displayed here in the plot insets, (a) CO/ = — gov2, (b) CO/ = 0, and (c) CO; = go\/2. The additional red dashed line plots 
correspond to a QD with no non-radiative dephasing (j p /2k = GHz). 



(the leftmost and rightmost peaks in the reflection spectrum) 
and for Yp=0. Higher levels of non-radiative dephasing lead 
to a washing out of the antibunching, though the field exhibits 
sub-Poissonian statistics on a time scale of ~ 1/ (K7- +y±)/2. 
For Cu/=0 (Fig. Qjjb)), the calculation predicts photon bunch- 
ing. This occurs because the resonance at zero detuning only 
appears when at least one photon is in the cavity and is cou- 
pled to the QD, which then allows for an additional photon 
at this frequency to be stored in the cavity and reflected. The 

(2) 

high-frequency oscillations in gccw,ccw(t), an d resulting nar- 

(21 

row super-Poissonian central peak about g C cw,ccw(0), are a re- 
sult of interference effects created by beating between the two 
cavity modes that are excited in this case (note that the central 
resonance peak in the cavity-QD spectrum is predominantly 
atomic-like, and excitation through the optical channel effec- 
tively excites the two detuned peaks, which are primarily pho- 
tonic in nature). 

Finally, we consider the parameter set explored in Fig. [TT] 
where p=0 so that only the QD couples the clockwise and 
counterclockwise modes together. For the Rabi-split peaks 
centered at CO/ = ±gov / 2, we see strong antibunching, as to be 
expected for a single QD coupled to a single cavity mode ex- 
cited on resonance with the Rabi-split peaks. At CO/ = 0, there 
are just minor oscillations about g^ (x)=l due to the weak and 
transient coupling of the resonant cavity mode with the QD. 



Comparison of this example with that of Fig. [T3l illustrates 
well the added system information gained by studying inten- 
sity correlations of the scattered light. Although both systems 
look very similar when studying the amplitude of light trans- 
mission and reflection intensity under weak-driving, the inten- 
sity correlations provide information about the spatial position 
of the QD relative to each of the standing wave cavity modes 
and the relative strength of go to p. 

Up to this point we have considered only the two-time cor- 
relation function for the counterclockwise propagating cav- 
ity mode. Correlation functions for the clockwise propagat- 
ing mode and standing wave cavity modes can be determined 
through formulas analogous to eqn. ( |23l , and can provide 
further insight into the appropriateness of the standing wave 
mode picture. Figure Q3] shows the results of two-time inten- 
sity correlation calculations for the set of parameters consid- 
ered in Fig. [14] where now we have plotted the intensity auto- 
correlation function for the standing wave modes. The results 
are consistent with the standing wave mode picture of atom- 
cavity interaction: mode swi is spatially aligned with the QD, 

(21 

and hence Wj (x) shows significant photon antibunching 
at the Rabi-split frequencies CO/ = ±goy/2 and off-resonance 
bunching[Q~2| at CO/ = (top plots in Fig. H3a)-(c)), while 
standing wave mode SW2 is spatially misaligned from the QD 

(21 

and gsw 2 ,sw2(t) i s essentially unity for all drive frequencies 
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FIG. 15: Normalized second order correlation functions (solid blue line) gi w \ , SWi (x) (top plots) and g sv j 2 ^ (t) (bottom plots) for the parameters 
{go i P, Kj- , K e , Y|| , y p }/2iz = {6,0, 1.2,0.44,0.16,0.7} GHz, Aco ac = 0, and whose transmission/reflection spectrum, originally shown in Fig. II II 

is re-displayed here in the plot insets, (a) CO; = — go\/2, (b) CO/ = 0, and (c) CO/ = go\/2. The additional red dashed line plots correspond to a 
QD with no non-radiative dephasing (y p /2n = GHz). 



(bottom plots in Fig. [T5la)-(c)). 

In addition to the autocorrelation calculations presented 
thus far, there are a number of other investigations of non- 
classical behavior within this system that may be of interest. 
For example, mixed-mode correlation functions can give in- 
sight into entanglement between the two cavity modes and 
the potential for generating non-classical states such as those 
employed in studies of the Einstein-Podolsky-Rosen paradox 
11551 15711 . Squeezing, which has been studied in the con- 
text of the Jaynes-Cummings system by a number of authors 
JM EE El, 113, EH, is also a potential topic for further study. 
The strong coupling of two cavity modes to a single QD, in the 
presence of background-mediated intermodal coupling, may 
yield important differences from previously studied systems. 
Furthermore, generating squeezed light or other non-classical 
fields in a microchip-based geometry could be of technolog- 
ical importance. In the Appendix, we present some prelimi- 
nary calculations on the above topics which indicate the de- 
gree to which such non-classical behavior may be exhibited 
in these devices. A more authoritative treatment of these 
topics requires a systematic investigation of different param- 
eter regimes for {go>^,^T,T^e,J\\,Jp}, driving field strength 
and frequency, and excitation channel (coupling to the cav- 
ity mode versus coupling to the QD directly), and is beyond 
the scope of this paper. 

VI. SUMMARY 

We have extended the standard quantum master equation 
model for a two-level system coupled to the mode of an elec- 
tromagnetic cavity to better reflect the situation that occurs in 
realistic semiconductor microdisk cavities. In this model the 
quantum dot, still treated as a two-level system, is coupled 
to two cavity modes corresponding to clockwise and counter- 
clockwise propagating whispering-gallery modes of the disk. 



These two modes are in turn passively coupled to each other 
through surface roughness, characterized by a backscatter pa- 
rameter p. We examine the steady state behavior of the system 
for differing regimes of P, the QD-cavity coupling rate go, the 
cavity decay rate Kt, and the quantum dot dephasing rate y_|_. 
In particular, we consider conditions for which standing wave 
cavity modes form, how the magnitude of the different system 
rates and the phase of p determine the nodes and antinodes 
of the cavity modes with respect to the quantum dot, and the 
resulting QD-cavity coupling. It is anticipated that this anal- 
ysis will be useful in the interpretation of experimental spec- 
tra from a waveguide-coupled whispering-gallery-mode mi- 
crocavity strongly coupled to a single two-level system such 
as the exciton state of a self-assembled quantum dot. 
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APPENDIX: NON-CLASSICAL CORRELATIONS AND 
SQUEEZING 

1. Mixed-mode correlation function 

We reconsider the parameter set of Fig. [14] where p=0, 
so that only the QD is coupling the two cavity modes, and 
calculate the mixed-mode correlation functions g C w,ccw{t) and 

(2) 

gswi,sw2\v f° r tne same set of parameters. We focus on solu- 
tions for t=0, which can be obtained entirely from the steady- 

(2) 

state density matrix, and examine the behavior of g a b {Q) as a 
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FIG. 16: Second order correlation functions for the parameters {go, P,Kj-, K e ,y|| , y p }/27t = {6,0,1.2,0.44,0.16,0.7} GHz, Aco fle = 0, and 
whose transmission/reflection spectrum, originally shown in Fig. [TT] is re-displayed here in the plot insets, (a) (gew,ccw(0)) 2 (solid blue line), 

gcw,cw(0) *gccw,ccw(0) (dashed red line), (giw ljJW2 (0)) 2 (dash-dotted green line), and g^ ltSWl (0) *giwl,sw 2 (0) (dotted black line) as a function 

(2) (2) (2) (2) (2) (2) 

of laser-cavity detuning, (b) g c J,ccw{t), gcw,cw{t), gccw,ccw(t), gsWi,sw 2 W, giwuswAt), and gm 2 ,sw 2 {t) ror excitation at C0//2JI = -6 GHz, and 

when both cavity modes are initially in the vacuum state and the QD is in its ground state. 



function of driving frequency. The Cauchy-Schwarz inequal- 
ity 



(*$f<&*$ (A-l) 

is violated when non-classical correlations exist between the 
two modes a and b lf55l l56ll . As we see in Fig. [ToTa). this 
inequality is violated for the traveling wave modes at par- 
ticular choices of CO/, so that non-classical correlations be- 
tween the two modes can occur in this system. On the other 
hand, for the standing wave modes, no quantum correla- 
tions exist, as mode sw2 is not coupled to the QD, so that 

(2) (2) 

gsJ u sw 2 (0) =8sJ 2 ,sw 2 (0) = 1 for all CO/. The transient (i.e., 
non-steady-state) behavior of the mixed and single mode cor- 
relation functions are shown in Fig. [ToT b). Here we plot 

8cw,ccw(t\ 8swi,sw 2 \^)> 8cw,cw(t\ 8 ecw,ccw (0 s 8swi ,swi and 

8sw 2 ,sw 2 (t) for CO//27I = —6 GHz (where the Cauchy-Schwarz 
inequality is nearly maximally violated), with the t = ini- 
tial state consisting of both cavity modes in the vacuum state 
and the QD in its ground state. The calculations indicate that 
steady-state behavior is achieved after ~ 1 ns, corresponding 
to the system's average decay time (~ 1 / (kt +J±_)/2), with 
violations of the Cauchy-Schwarz inequality for the traveling 
wave modes occurring after only ~0.2 ns. 

2. Squeezing 

As has been observed by several other authors in studies 
of single mode cavity QED 05 8l 16011 . squeezing in the field 
quadratures can occur, however, with the amount of squeez- 
ing typically small (< 20%) unless large intracavity photon 
numbers (> 10) are achieved II 6411 . The basic reason for this is 
that the nonlinear interaction that generates squeezing in the 
Jaynes-Cummings system is that of the electromagnetic field 



coupling to a saturable oscillator (the QD); this implies that 
the intracavity field has to be strong enough for QD saturation 
effects to be appreciable. For our system, situations where 
only one of the standing wave modes is coupled to the QD 
(as in Figs. [T2]and[T4] for example) essentially reduce to that 
of the single mode cavity QED case, and we expect qualita- 
tively similar behavior. A perhaps more interesting example 
to study is that of Fig. [13] where both standing wave modes 
are equally coupled to a QD. To achieve a reasonable intra- 
cavity photon number, we increase the input driving field by 
approximately three orders of magnitude over the weak drive 
fields we have used up to this point, to a level of ~ 30 pho- 
tons/ns, so that the average intracavity photon number (Fig. 
fTTl b)) peaks at a value of ~ 1 . This results in the transmission 
and reflection spectra shown at the top of Fig. [TTl' a). In com- 
parison to the transmission/reflection spectra calculated in the 
weak driving limit in Fig. |8la), we now begin to see asym- 
metries in the transmission dips (reflection peaks) that are as- 
sociated with multi-photon transitions to excited states in the 
Jaynes-Cummings spectrum and QD saturation effects. In this 
calculation, we are unable to numerically study higher drive 
strengths due to the resulting large system size for these two- 
mode cavities. To access higher driving fields using the same 
computational resources, adopting a wavefunction-based ap- 
proach (i.e., the quantum Monte Carlo method) is one possi- 
bility. 

We next consider fluctuations in the steady-state, internal 
cavity field (squeezing in the external field, which can be in- 
vestigated through the spectrum of squeezing [62], for exam- 
ple, are not considered here but may be of future interest). 
First, we look at fluctuations in the photon number in mode i 
by calculating the Mandel Q parameter [51] 



n( *s Var(a[fl/) - (aja/) 

2K) = TTfTT 

(«/«;) 
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FIG. 17: Non-classical properties of the microdisk cavity-QD system under more intense driving (~30 photons/ns drive power), and with 
system parameters {gOj^KT-jK^Yi^Ypl^Jl = {6, i9. 6, 1.2,0.44,0.16,0.7} GHz. (a) Cavity transmisson/refiection as a function of driving 
field frequency CO/, (b) total (solid black line), CW mode (dashed blue line), and CCW mode (dot-dashed red line) intracavity photon number 
as a function of CO/, (c) Q parameter for the intracavity photon number in the (top) traveling wave WGMs and (bottom) standing wave WGMs. 
(d) Q parameter for the X 1 and X 2 quadratures of the (top) cw and (bottom) ccw traveling wave mode. 



X 



X 



(A3) 



where for some operator O, Var(O) = ((9 2 ) — {O) 2 . Fig- tures. For mode i, we define the quadrature operators x/' 2 by: 
ure \T7\ c) shows the calculated Q parameter as a function of 
driving field frequency for the cw/ccw traveling wave modes 
(top) and swi/s\V2 standing wave modes (bottom). These 
plots show Q < for certain driving frequencies, indicating 
that sub-Poissonian photon number statistics can be achieved, 
though the level of non-classicality is small (~ 5%). A calcu- 
lation of <2(«cw) (f ) for <M/=0 ( with both cavity modes initially 
in the vacuum state and the QD in its ground state) indicates 
that slightly higher levels of non-classicality (Q ~ —0.1) can 
be achieved before steady-state is reached. Additional pre- 
liminary calculations using a quantum Monte Carlo method to 
access higher drive strengths have been performed, and show 
that Q can continue to decrease for larger driving fields. For 
drive strengths of ~ 300 photons/ns (corresponding to an av- 
erage total intracavity photon number ~1 at C0/=0), 2(«cw) 
can reach —0.35 in its transient (non-steady-state) behavior. 
Similarly, one can examine fluctuations in the field quadra- 



The corresponding Q parameter for quadrature j of mode i is 
then JUl: 



Ql 



Var(l/)-0.25 



0.25 



(A.4) 



From Fig. [T7T d). we see that small amounts of quadrature 
squeezing in the cw/ccw modes are apparent for the condi- 
tions considered. 
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